Coherence scale of the Kondo lattice 
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It is shown that the large-N approach yields two energy scales for the Kondo lattice model. The 
single-impurity Kondo temperature, Tk, signals the onset of local singlet formation, while Fermi 
liquid coherence sets in only below a lower scale, T*. At low conduction electron density n c ("ex- 
haustion" limit) , the ratio T* /Tk is much smaller than unity, and is shown to depend only on 
n c and not on the Kondo coupling. The physical meaning of these two scales is demonstrated by 
computing several quantities as a function of n c and temperature. 



How coherent quasi-particles form in the Kondo lattice 
has been a long-standing issue. For a single impurity, 
there is a single scale Tk below which the local moment 
is screened and a local Fermi liquid picture applies. Tk 
can be defined e.g. as the scale at which the effective 
Kondo coupling becomes large. All physical quantities 
(e.g. specific heat, susceptibility) obey scaling properties 
as a function of T/Tk in the limit where both T and 
Tk are smaller than the high-energy cutoff (bandwidth, 
denoted 2D in the following). In contrast, for a Kondo 
lattice, one may suspect that the physics is no longer 
governed by a single scale. Indeed, while magnetic mo- 
ments can still be screened locally for temperatures lower 
than the single- impurity Kondo scale Tk, the formation 
of a Fermi-liquid regime (with coherent quasi-particles 
and a "large" Fermi surface encompassing both conduc- 
tion electrons and the localized spins) is a global phe- 
nomenon requiring coherence over the whole lattice. If 
at all possible, it could be associated with a much lower 
coherence temperature T*. The situation is reminiscent 
of strong-coupling superconductivity in which local pair 
formation may occur at a much higher scale than T c , the 
scale at which long-range order sets in. 

As was originally pointed out by Nozieres JO , this issue 
becomes especially relevant when few conduction elec- 
trons are available to screen the local spins, i.e. in the 
limit of low concentration (n c « 1). In this "exhaus- 
tion" regime, two possibilities arise: i) either magnetic 
ordering wins over Kondo screening or ii) a paramag- 
netic Fermi liquid state still manages to form, but with 
a much suppressed coherence scale T* <C Tk- Nozieres 
has suggested in Ref. that T* ~ n c D for strong cou- 
pling J K /D > 1 (where T K ~ Jk), while T* ~ T K /D 
for weak coupling Jk <C D. Recently, several studies 
HHHHEEH nave addressed this issue using dynamical 
mean- field theory (DMFT) j|. The conclusion was that 
two scales are indeed present in the Kondo lattice, with 
the coherence scale with T* <C Tk in the "exhaustion" 
regime. Because the DMFT equations require a numer- 
ical treatment, no detailed analytical insight into these 



two scales was obtained, even though the validity of the 
estimate || T* cx T K /D was questioned M. 

In this letter, we solve this problem using the slave- 
boson approach, in the form of a controlled large-N so- 
lution of the SU(N) Kondo lattice model. This ap- 
proach has been extensively used in the past twenty years 
0,1^1111111]. Surprisingly, the issue of the coher- 
ence scale and the temperature dependence of physical 
quantities has not been discussed in detail in the exhaus- 
tion limit n c <C 1 (See, however, Ref. jl6]]). We find 
that the large-N approach provides a remarkably simple 
and physically transparent realization of the two-scale 
scenario described above. Furthermore, because of its 
simplicity, it allows for an analytical calculation of the 
coherence scale, which is found to disagree (for weak cou- 
pling) with Nozieres' estimate in Ref. ^| (while it agrees 
with it at strong coupling). We also calculate the tem- 
perature dependence of several physical quantities and 
find remarkable agreement with the more sophisticated 
(and demanding) DMFT approach. 

We consider the Kondo lattice model (KLM): 
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In this expression, the spin symmetry has been extended 
to SU (N) (<t = 1, • ■ • , N) and the local spins will be con- 
sidered in the fermionic representation: Sf a = f^fia' — 
5 aa >/2, with the constraint J2a ffofia = N/2, Standard 
methods Jl(| are used to solve this model at large N: a 
boson field B^t) (conjugate to the amplitude f^Cia) 
is introduced in order to decouple the Kondo interaction, 
and the constraint is implemented through a Lagrange 
multiplier field A^(r). For N = oo, a saddle point is 
found at which the Bose field condenses (Bi(r)) = r 
and the Lagrange multiplier takes a uniform, static value 
(iAi(r)) = A. Two quasiparticle bands ui^ are formed, 
solutions of (uo + X)(u> + (i — £fe) — r 2 = 0. Changing vari- 
ables to wjr, the three saddle-point equations can be cast 
in the compact form: 
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1 1 n c 
Jk 2' 2 



dLonp(uj)po(u) + p 



uj + A' 



F*(n c ) = exp 



p (£F +w) - Po(cf) 



cj + A' (w + A) 2 ' 



(2) 



Here, Uf is the Fermi function, n c /2 is the conduction 
electron density per spin colour, and po(e) = J2k ^( e — e fc) 
is the non-interacting density of states. 

In the large-N approach, the onset of Kondo screening 
is signaled by a phase transition at a critical temperature 
T/f, below which r(T) becomes non-zero. The equation 
for T K is 2/J K = J de Po {e) tanh[(e - p )/2T K ]/(e - Mo), 
with po the non-interacting chemical potential (at T = 
Tk)- This equation coincides with that for the single- 
impurity case: Kondo screening of individual local mo- 
ments starts taking place in this approach precisely at the 
single-impurity Kondo scale. We have derived an explicit 
expression for Tk in the weak-coupling regime Jk *C D: 



T k = De-V^W ^1 - ( Cp / D) 2 F K (n B ) , (3) 
' f D ~ CF duj p (e F + (j) - Po(eF) 

'-(£)+ep) 



F K (n c ) = exp 



2po(£f) 



I, (4) 



where £f is the non-interacting Fermi level (given by 
n c /2 = depo(e))- This expression is valid for an even 
d.o.s po(— e) = jOo(e) which vanishes outside the interval 
— Z? < e < +D. The factor Fr-, equal to unity for a con- 
stant density of states, varies smoothly with n c (or cf)- 
In contrast, the prefactor [1 — (e F / D) 2 ] 1 / 2 vanishes in the 
low-density limit n c — > 0, and suppresses Tk. 

We now consider the low-temperature limit T <C Tr-, 
in which the large-N approach leads to an extremely 
simple Fermi liquid picture fTcj ] . It is somewhat over- 
simplified in that the conduction electron self-energy 
Ti c {k,uj) — r 2 /(uo + A) is purely real (finite lifetime ef- 
fects are absent at N = oo) and momentum indepen- 
dent. Even so, it captures some of the most impor- 
tant features of the problem. The zero-frequency shift 
S c (w = 0,T = 0) = r(T = 0) 2 /A(T = 0) is precisely 
such that the Fermi surface has a large volume contain- 
ing both conduction electrons and local spins. Indeed, 
adding the last two saddle point equations in (^|), one 
obtains: p — £ c (w = 0, T = 0) = e F where e F is the 
non-interacting value of the Fermi-level corresponding to 
(n c + l)/2 fermions per spin color. As detailed below, all 
physical quantities at T = are directly related to a sin- 
gle energy scale, proportional to the boson condensation 
amplitude T* — r 2 (T — 0)/D. It is possible to derive an 
analytical expression for T* in the weak coupling limit, 
which is valid for a general (bounded) density of states 
and arbitrary density n c . This expression, which has ap- 
parently not been reported before, reads: 

T* = De- 1/jKpo(eF) (1 + e F /D) ^-F*{n c ), (5) 



Po(£f) 



, (6) 



where Aep = ep — ep- As Fk, F* varies smoothly with 
n c |l7]]. The total density of states at the Fermi level 
p(u = 0) = p cc {0) + p ff (0) is given by: p(0) = p (e F )/Z c , 
with Z c the conduction electron quasiparticle residue 

l/Z c = 1 - d£ c ( W )/cM w=0 = 1 + ( £ f - ep) 2 Ao- In 
the weak coupling limit, all physical quantities at T = 
are directly related to p(0) ~ p (e^)(Ae F ) 2 /(T*D), and 
are thus renormalized by the ratio T*/D. For instance, 
the /-electron susceptibility and the specific heat coef- 
ficient (per spin color), are given by \f K p(0) an d 
7 = 7r 2 p(0)/3. The Drude weight may also be computed 
with the result D R oc r*Dp (eF)/(Ae F ) 2 - 

The physical content of these expressions and of 
Eqs.(|[^) is that two different energy scales are rele- 
vant for the Kondo lattice model: one (Tk) is associ- 
ated with the onset of local Kondo screening; the other 
(T*) with Fermi liquid coherence and the behaviour of 
physical quantities at T = 0. These two scales have the 
same exponential dependence on Jk/D at weak coupling, 
but very different dependence on the conduction elec- 
tron density in the "exhaustion" limit n c C 1, in which 
T* <C Tk ■ This is in qualitative agreement with Nozieres 
proposal in [|| , but not with his estimate T* /Tk oc Tk / D 
(which would thus depend on the coupling). We find 
that the ratio T* /Tk depends only on n c in this limit, 
in a manner which reflects the behavior of the bare d.o.s 
Po(e) at the bottom of the band. For a d.o.s vanishing as 
(e + D) a , Eqs.§|) yield T*/T K oc nV {2[1+a)] . 

Fig. |] displays the rt c -dependence of Tk and of the in- 
verse of the /-electron susceptibility, obtained by solving 
the large-N equations, both for a single-impurity and for 
the lattice. It is seen that Tk and l/x/(T = 0) are iden- 
tical energy scales for the single impurity case, but have 
very different density dependence for the lattice model 

1 /3 

(T* /Tk oc n c at low density for the semi-circular d.o.s 
used here). Near n c = 1, Tk falls below 1/x reflecting 
the vanishing \ for the Kondo insulator. These curves 
are remarkably similar to those obtained by Tahvildar- 
Zadeh et al. in their Quantum Monte Carlo studies of 
the periodic Anderson model in infinite dimensions || 
(see also ||- [||). In the inset of Fig. [j], we plot the di- 
mensionless number TkXi{^) as a function of IuTk/D. 
This number goes to a universal value at weak coupling 
in the single-impurity case, while it has intrinsic density 
dependence for the lattice. This shows that Tk is not the 
appropriate low-temperature scale, especially for n c — » 0. 
The inset in Fig. || shows the effective mass ratio m/m* 
(computed from the specific heat) and the Drude weight, 
as a function of n c . This vanishes in the limits n c — > 
and n c — > 1. The effective mass m* is proportional to 
e +i/ JKpo(e F ) ^ with a density dependent prefactor which 
diverges as n c — > and vanishes at half filling. 
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FIG. 1. Solid line: density-dependence of Tk (denned 
as the slave boson condensation temperature). Squares: 
l/4x/(T=0) for a single impurity. Dashed line: l/4x/(T=0) 
for the lattice. Jk / D=0.75. Inset: 2V-dependence of 
TxX/(0) for the impurity (left scale, solid symbols) and the 
lattice (right scale, open symbols) for n c =0.1, 1.0 (impurity) 
and n c =0.1, 0.5, 0.9 (lattice) from top to bottom. All plots 
are for a semi-circular d.o.s. 



We have also studied the behavior of several physical 
quantities as a function of temperature, by solving nu- 
merically Eqs.(|2|). Fig.|| shows the product Tkx{T) for 
the Kondo lattice, for several values of function 
of T/T K . All the curves merge at T/T K = 1, where the 
boson decondenses and Xf reaches the free spin value 
( Xf (T) = 1/4T iorT >T K ). No universal scaling func- 
tion of T/Tk describes the temperature dependence in 
Fig. ||, in contrast to the single-impurity case. Plotting 
\(T)/x{T — 0) does not restore scaling, since qualitative 
differences in the T-dependence are seen for different den- 
sities. Fig. |](a) shows the T-dependence of the entropy. 
A linear behavior S(T) = 7T is found at low tempera- 
ture for all densities n c ^ 1. The slope 7 oc p(0) decreases 
with increasing density as does the temperature scale Tp 
up to which SiT) is linear. Tp is of the order of T* 
(<C Tk) at low n c , while it can be estimated by compar- 
ing 7 T to e~ TK / T as T* ~T K /\ ln(l - n c )| < T K T* 
for n c — ► 1. For n c c± 1, Tp is a better estimate of the 
coherence scale than T* itself. 

At low densities, after a steep initial rise, S(T) re- 
mains close to ln2 up to Tk- This can be interpreted 
in terms of the strong-coupling picture discussed below. 
At higher densities, most of the variation takes place 
in the range T* < T < T K . The specific heat C(T), 
shown in Fig.g(b), has a two-peak structure |L8|: the 
peak at Tk signals the onset of Kondo screening and 
appears in this mean-field description as a discontinu- 
ity of C(T). The second peak, at Tp, signals the Fermi 
liquid heavy-fermion regime. As n c increases, weight is 
gradually transferred to the high-temperature peak until 
the low temperature peak disappears completely in the 
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FIG. 2. T-dependence of the /-electron susceptibility for 
the lattice for n c =0.01, 0.2, 0.5, 0.8 and 0.9, from top to 
bottom. Inset: Inverse effective-mass(left scale, full line) and 
Drude factor (right scale, dashed line) as functions of density. 
In all plots, J K /D = 0.75 
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FIG. 3. Entropy (a) and specific heat (b) for the lattice 
model for several values of the density. Jk I D=0.75. 



Kondo insulator limit. 

In the strong coupling limit Jk/D ~S> 1, the large-N 
results support the physical picture proposed by Lacroix 
p9| and further discussed in In this picture, the con- 
duction electrons bind to n c spins, forming singlets be- 
low T ~ Jk ~ Tk, the binding energy of a singlet. The 
remaining 1 — n c "bachelor spins" behave as itinerant 
fermions subject to a constraint of no double occupancy. 
The hopping integral of the resulting effective infinite-U 
Hubbard model is t c ff = —t/2. The (hole-like) sign of this 
matrix element implies that these 1 — n c fermions have 
a Fermi surface volume corresponding precisely to n c + 1 
particles. In the exhaustion limit n c <C 1, one has effec- 
tively a weakly doped U — 00 Hubbard model. Solving 
Eqs.(||) at strong-coupling yields a quasiparticle residue 
Z c oc n c , and hence a coherence scale T* ~ n c D cor- 
responding to a Brinkman-Rice estimate for this doped 
Mott insulator 101. Notice, however, that at finite N 
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FIG. 4. n c (e) for T/T K =1.Q (solid line), 0.5 (dashed line), 
0.005 (dotted line).ro c =0.15 Inset: Spectral density T(u>) for 
TjT K =\ (solid line), 0.5 (dashed line) and 0.25 (dotted line). 
For T < Tk there is a 5-function peak in the gap not shown 
for clarity. n c =0.2. In all plots Jk / D=0.75. 



this uniform solution may become unstable to magnetism 
or phase separation of localized singlets and unscreened 
spins. This picture sheds some light on the T-dependence 
of the entropy and specific heat in the exhaustion limit 
reported above. As the system goes through Tk, it looses 
the magnetic entropy (~ n c In 2) of the n c bound spins. 
The remaining entropy (~ (1 — n c )ln2), is lost below 
T* . The two peaks of unequal weight in the specific heat 
reflect these processe s pc| . 

We display in Fig.|]the distribution function of the 
conduction electrons n c (ek) = (ctck). Very close to 
Tk, n>c(tk) has the shape of a Fermi function centered 
around the non-interacting Fermi level ep, with a small 
thermal broadening of order Tk (<SC D in weak cou- 
pling). As T is reduced below Tk, weight is transferred 
to scales of order e^, the interacting Fermi level asso- 
ciated with the large Fermi surface. The feature at ep 
broadens as T is decreased. When Fermi liquid coher- 
ence establishes at T ~ T F , a discontinuity (of ampli- 
tude Z c ) develops at e = e^. Finally we note that, in the 
large-N approach, the Kondo lattice model can be exactly 
mapped onto an effective single-impurity model coupled 
to a self-consistent bath of electrons. This mapping holds 
more generally for any approach in which the conduc- 
tion electron self-energy depends only on frequency, such 
as dynamical mean-field theory ||. For a semi-circular 
d.o.s, the Green's function of the self-consistent bath is 
g (icu n ) = [iuj n +n-D 2 G c (iuj n )m~ l . The inset of Fig. | 
displays the continuous evolution of the spectral density 
T(u) = — ImQo(u + i0 + )/ir as T is reduced below Tk- 
Above Tk, it coincides with the non- interacting d.o.s po, 
while the bath is split into two bands below Tk- The 
width of the lower band depends on n c and is small for 
« c « 1, leading to the low coherence scale. The existence 
of a sharp gap separating the two bands is an artifact of 



the large-N limit (except at n c = 1). In more realistic 
treatments it is replaced by a pseudo-gap [||J^^1. 

In conclusion, we have reconsidered the time-honored 
large-N approach to the Kondo lattice model, with spe- 
cial emphasis on the "exhaustion" limit of low electron 
density. We showed that two energy scales appear for 
which we have obtained explicit analytic expressions : 
the Kondo scale associated with the onset of local Kondo 
screening and a much lower scale associated with Fermi 
liquid coherence. Physical quantities reflect the crucial 
role played by these two scales. In this approach, mag- 
netic ordering is suppressed by quantum fluctuations : 
more elaborate treatments such as DMFT must be used 
to assess whether the coherence scale can actually be 
reached or magnetic ordering sets in at a higher tem- 
perature. 
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